Kertesz on Fat Graphs? 
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Abstract 

The identification of phase transition points, /3c, with the percolation thresholds 
of suitably defined clusters of spins has proved immensely fruitful in many areas 
of statistical mechanics. Some time ago Kertesz suggested that such percolation 
thresholds for models defined in field might also have measurable physical conse- 
quences for regions of the phase diagram below /?c, giving rise to a "Kertesz line" 
running between Pc and the bond percolation threshold, Pp, in the M, f3 plane. 

Although no thermodynamic singularities were associated with this line it could 
still be divined by looking for a change in the behaviour of high-field series for 
quantities such as the free energy or magnetisation. Adler and Stauffer did precisely 
this with some pre-existing series for the regular square lattice and simple cubic 
lattice Ising models and did, indeed, find evidence for such a change in high-field 
series around f3p. Since there is a general dearth of high-field series there has been 
no other work along these lines. 

In this paper we use the solution of the Ising model in field on planar random 
graphs by Boulatov and Kazakov to carry out a similar exercise for the Ising model 
on random graphs (i.e. coupled to 2D quantum gravity). We generate a high-field 
series for the Ising model on <I>^ random graphs and examine its behaviour for 
evidence of a Kertesz line. 



1 Introduction 



The question of how to give a geometrical, percolation-hke interpretation of the thermal 
transition in the Ising model was finally resolved by Coniglio and Klein [jl|, although the 
idea is already implicit in the work of Fortuin and Kasteleyn P| on a correlated bond- 
percolation model for the g-state Potts model and in many ways is a realisation of Fisher's 
earlier ideas on critical droplets 0. The key insight was to construct spin clusters in which 
like spins were joined with probability p = 1 — exp(— 2/3), where f3 = l/fc^T is the inverse 
temperature. These then percolated at the correct critical temperature and gave the 
correct thermal exponents. The definition of such stochastic clusters has led directly to 
the development of cluster algorithms and similar geometrical representations in many 
statistical mechanical models. 

Since Coniglio-Klein clusters are intended to access the thermal phase transition of 
the Ising model they are not usually formulated in the presence of an external field, since 
in this case the free energy becomes analytic and the thermal transition is washed out. 
However it is clear that one can still build such clusters when the external field is non- 
zero. For a given external field value H these clusters will percolate at some temperature 
Tx{H), and varying H then traces out the so-called Kertesz line shown in Fig. 1 in the 
M,u plane, where M is the magnetisation and u = exp(— 4/5). For if = we have 
the usual thermal Ising transition at Uc, whereas for if — ^ oo we recover standard bond 
percolation for the lattice in question at Up. 

Although a percolative transition exists across the Kertesz line no singularities are ex- 
pected in the thermodynamic quantities related to the Ising model, since these can only 
appear for H = 0. Nonetheless, something physical is going on since the percolative tran- 
sition is signalled by "wrong-sign", say down, spin clusters losing their surface tension. 
This begs the question of whether the Kertesz line leaves any discernible signal in mea- 
surable quantities. This issue was discussed by Kertesz himself [Q, and later investigated 
numerically by Stauffer and Adler |^. If one considers the number of down spin clusters 
ns{u, H) in a non-zero up field H then standard nucleation theory leads to 

logns{u,H) = slogy - sT (1) 

below the Kertesz line, where y = exp{—2H) and a = 1 — 1/d. The second term arises 
from the surface tension F of the droplets and is a direct extrapolation of the H = 
{y = 1) result. On the other hand above the Kertesz line a zero surface tension expression 
of the form 

\ogns{u, H) = slogy + s\og\{u, H) (2) 
where A < 1 is postulated. This means that cluster numbers will decay as 

n, = exp{-Ts^~^/'^)y' (3) 

below the Kertesz line whereas they will decay as y^A^ above the line. 

Physical quantities such as the magnetisation can be expressed in terms of n^, 

M = l-2Y,sns{u,H), (4) 
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so we can see that they can be expressed in terms of a (large-field) expansion in y using 
equs. (|I]) and (0). Since equ. (0) shows that y is effectively replaced by yX above the 
Kertesz line, one might expect the series for M (and others) to converge up to some y > 1 
(since A < 1). 



2 Square Lattice Reprise 

In order to investigate the possible consequences of a Kertesz line Adler and Stauffer 
analysed the radius of convergence of some long-standing high-field series for the free 
energy and the magnetisation of the Ising model using Pade approximants 0. We briefly 
review their results here for comparison with the random graphs in the next section. 
Taking the series for the free energy F or its derivative, we have 

F = J2Ls{u)y% 

s 

where the so-called high-fleld polynomials Ls{u) were calculated up to order 15 for the 
square lattice in 

The series for F, g and the magnetisation M were subjected to a standard unbi- 
ased Dlog-Pade analysis in order to estimate the radius of convergence. If we assume a 
singularity of the form F ~ (y — yc)~^, then 

^^ogF = ^ (6) 

dy y-Vc 

and an [L/M] Fade approximant to ^logF (or or M), 

^ ^ ^ Quiy) qo + qiy + ---qMy''' ^ ' 

would be expected to show the pole at yc as a zero of Qniy)- For the square lattice for 
instance Ur = 3 - 3^2 = 0.17157 . . . and Mp = 1/4, so we would expect to see a change in 
yc at, or around, Up. 

We show the results of such an analysis on ^ for the estimated radius of convergence 
as a function of u in Fig. 2. We can see that there is evidence for a change in the 
behaviour of yc around m = 1/4. In the plotted data jumps are present in the [7/7] and 
[6/7] approximants at Wp, and the [5/5] and [6/6] approximants cease to give a real pole 
term for some u ~ Up. In Fig. 2 we have plotted the [5/5] and [6/6] approximants up to 
the largest real pole value obtained. As shown in [Q the behaviour for other approximants 
and for F and M is similar - in all cases there is evidence of an increase in the estimated 
yc around m = 1/4 or the approximant's estimates of yc become complex. 

One detail that should be noted is that both here and in the random graph case inves- 
tigated below there is a tendency for the approximants to throw up spurious cancelling 



2 



pole/zero pairs, so one must be careful to look at both the numerator and denominator 
of the approximant to make sure that one is determining the "real" pole term. This also 
introduces a certain degree of subjectivity into proceedings since the poles and zeroes are 
only equal to within numerical accuracy. All the series calculations were done with exact 
arithmetic, with fixed (but large) precision only being resorted to in the final stage of 
obtaining the roots of the Pade numerators and denominators. 

An analysis of the series for the simple cubic lattice Ising model, available up to 
order 13, in |^ also gave very similar results, with signs of an increase in the estimated 
radius of convergence, or instability in the approximants, in the vicinity of the percolation 
temperature for the lattice. The usual caveats of course apply to any such discussion, since 
the series involved are quite short (and of some considerable vintage), but the cumulative 
evidence of the various approximants and the results on different lattices was felt in |^ 
to lend support to the Kertesz scenario. 



3 Fat Random Graphs 

Given the evidence for a Kertesz line in the Ising model on the square and cubic lattices it 
is tempting to look elsewhere for the phenomenon. In order to carry out such an analysis 
for other models (or lattices) we require a high-field expansion, which are in general rather 
hard to come by. In what follows we discuss obtaining such an expansion for an exact 
solution to the Ising model in field on planar (fat) random graphs which was derived by 
Boulatov and Kazakov |^, They considered the partition function for the Ising model 
on a single planar graph with n vertices 

^singie(G", /?,//) = ^ exp [(3Y.Gla,a, + HY,aA , (8) 

then summed it over some suitable class {G"} of n vertex graphs (e.g. or <l>^ random 
graphs) resulting in 

Zn = 51 Zsing\e{G"', P, H) , (9) 

before finally forming the grand-canonical sum over differently sized graphs 



n 



Zn, (10) 



where c = u^^"^ = exp(— 2/3). This last expression could be calculated exactly as matrix 
integral over N x N Hermitian matrices. 



Z = - log J V(pi V(p2 exp (^-Tr 



^(</.? + </.^)-c0i02-|(e^0t + e-^0^) 



where the N —>■ oo limit is to be taken to pick out the planar diagrams and we have used 
the potential appropriate for $^ (4-regular) random graphs. 
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When the integral is carried out the solution is given by 

„ 1 , z 1 dt , ^ 1 dt , 

Z = -\og / -g{t) + ^^ / -g{tf. 

2 g g Jo t 2g^ Jo t 

where the function g{z) is 

1 ^ 3z(yV2 + ^-1/2 



g(z) = 3C^Z-^ + Z 



3z) 



2) 



(1 - 9z 



2\2 



(12) 



(13) 



As we can see this is a solution in field for the Ising model on random graphs, since g{z) 
has been obtained in full generality for y = exp{—2H) ^ 1. 

In order to generate a high-field series from equ. (|l^) we revert the series for g{z) to 
get a series z{g). This can then be used in equ. ([T^ ) in order to obtain an expansion for 
the partition function in powers on g (i.e. in the number of vertices). As we noted in |^ 
it is only necessary to consider the first of the terms in Z = Zi + Z2 + with 



Zi 
Z2 
Zs 



2 9 
1 



dt 



g Jo t 



git), 



(14) 



1 



dt 



/ -9ity, 



2g^ Jo t 



since Z^ = J2n (^k^n{u)g^ where An{u) is identical for all the Zj^. 

This can be traced back to the generic expression for the partition function of any 
Hermitian matrix model, which is of the form 



z= /'rfe(i-Oiog(/(0) + ---, 

Jo 

where in the case of the Ising model on <l>^ graphs / is the solution to 



(15) 




2^Qgli_l^' + y-' 



-1/2 



c 




(16) 



The expression in equ. (|l^) is obtained from this by defining z = 2gf /c and integrating 
by parts. One finds that the coefficient of (7" in Zi should be {n + l){n + 2) times the 
full value obtained from expanding Zi + Z2 + Z3 from these considerations, which can 
be confirmed by comparison with the results presented for low orders in |T0, For 
example, the coefficient of g^ in |[T^ is found to beQ 



which is 



12 X Zi = 12 X {3/2)c-^y-^[9 + (16c^ + 2c'^)y + V 



;i7) 



^ Where g 
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as expected. The high-field polynomials come from an expansion of the free energy F 
rather than the partition function Z, so there is one further step. 

We choose to consider ^ for calculational convenience rather than directly taking 
the logarithm to obtain F, since the operation of taking a logarithm proved to be rather 
memory intensive. Since the free energy per site is defined as (absorbing a factor of —f3) 

F= lim - logZn, (19) 

n— »oo fl 

we can approximate this by 

d 1 dZn 

rNJ 

dy nZn dy 

for some sufficiently large n (in our case we take n = 32) 
in P was in any case similar to that of F and M. 



(20) 

. The behaviour observed for ^ 



4 The Percolation Threshold on <l> Graphs 

Kazakov has already calculated the threshold on graphs, so it is a simple matter to 
repeat the calculation to obtain the result he stated {pcr = 2/3), but did not derive, for $^ 
graphs in ^ . We have chosen to work with $^ graphs rather than $^ graphs here because 
the expression for g{z) turns out to be rather more inconvenient with <l>^ graphs and 
reverting this in the manner of the previous section to get Zi is much more cumbersome 
than it is for $^ graphs. 

The idea is to consider the percolation problem as the g — > 1 limit of a g-state Potts 
model, where the Potts partition function is 

^Potts = Eexp(-/37^). (21) 

with H = —J2{ij){^(Ti,aj — 1), and the spins cTj take on q values. This Potts partition 
function, just as in the Ising case, can be expressed as the matrix integral over N x N 
Hermitian matrices 0j 

F = ^hgll[V^,expi-S) (22) 

i 

where 

S=lt<Pl-~^EM,-lt<Pl (23) 

i=l i<j i=l 

and we have used c = l/(exp(/3) + g — 2) for the Potts coupling and temperature to 
distinguish them from the Ising case. The apparent factor of two difference between the 
q = 2 version of c and the Ising coupling in the previous section and the consequent factor 
of two in the temperature scales is accounted for by the use of a aiaj interaction in the 
Ising model and a Scr^^aj interaction in the Potts model. 
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This in turn may be recast into a matrix external field integral by introducing a further 



matrix integration in the dummy N x N Hermitian matrix X, 



X 



J 2)0 exp (hXcp 




(24) 
(25) 



where h"^ = c. In the limit g — > 1, the percolative probability p = 1 — exp(— /?) is related 
to these parameters by /i^ = c = 1/p — 1. The mean number of percolative clusters per 
unit volume f{p) may be calculated from the quantity C{g,p), which is given in terms of 
the Potts free energy by 

lim|^ = C(^?,p). (26) 

q^l Cjq 

Evaluating C{g,p) in the saddle-point approximation gives 



C{9,P) 



1 ^ 



^logA(x* 



(27) 



where the x* are the saddle point values of the X eigenvalues and F^{g) is the standard 
one-matrix model free energy. As in the model the log A(a;*) term is what counts 
for the percolative critical behaviour, and this gives 



logA(x*) 



p{u)p{v)dudv 



log 




+ UV + v^) 



2 p 



+ log |m — t;! 

(21 



where p is the eigenvalue density for the one-matrix $^ model which has support on 
[—a, a]. The model will be critical when the argument of the first logarithm is zero 
within the region of integration, which first occurs when = = = = 2/3, so 
Per = 2/3 as announced in [0]. Since Pcr = 2/3 = 1 — exp(— /3p), allowing for the factor 
of two between Ising and Potts conventions gives exp(— /3p) = exp(— 2/?p) = 1/3, so on 
the graphs Up = exp(— 4/5p) = 1/9 to be compared with the 1/4 of the regular square 
lattice. 



5 Pade Approximants for Graphs 

The procedure is identical to the square lattice investigation: here we take Z32, calculate 

^F ~ (29) 
dy 32 Z32 dy 

as our estimate of -^F and then take unbiased Dlog-Pade approximants to this expression 
in order to determine yc- Again, exact arithmetic is used until the final stage of extracting 
the zeroes and poles of the Pade approximant is reached. 
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Since our starting series is rather longer than the square lattice case, we calculate up to 
the [15/15] approximant. In Fig. 3 we show the [13/13], [14/14] and [15/15] approximants 
plotted against u. The near-diagonal approximants and lower approximants on the 
graphs which we also calculated all give very similar estimates for and we have not 
shown them for clarity. 

It is immediately obvious from Fig. 3 that although there is evidence for a change in 
the behaviour of the high-field series for u > Up, it is rather less clear cut than for the 
square lattice series. We have emphasised the fact by plotting the [7/7] square lattice ap- 
proximant over the $^ results to allow a direct comparison. Indeed, without the prejudice 
of a previous calculation of Up it would be fair to state that although there are some signs 
of an increase in the estimates of yc in the vicinity of = 1/9 the variation is strongly 
apparent only for larger u values. 

Our calculation of the series for $^ graphs is for a set which includes various degen- 
eracies such as self-energy and "setting-sun" diagrams in the ensemble of graphs, and the 



inclusion of such graphs has been found empirically to reduce finite-size effects |]12| by 
comparison with more restricted ensembles of graphs. It is of course possible that any 
ensemble of random graphs may be subject to greater finite-size corrections than a square 
lattice so it is conceivable that this might explain a weaker signal for the Kertesz line 
than on the square lattice. 



6 Conclusions 

We reviewed the existing results for the existence of the Kertesz line for the Ising model 
on the square lattice, which lend support to the suggestion that a change in the behaviour 
of the high-field series should be visible around the percolation threshold. We then made 
use of the Boulatov and Kazakov solution of the Ising model on planar random graphs 
to construct a high-field series for planar graphs and subjected this to the same Dlog- 
Pade analysis carried out on the (shorter) square lattice series. We found that there are 
clear signs in this series too for a change in behaviour, since the estimated values of the 
radius of convergence, yc, increase as m = exp(— 4/5) is increased, as seen in Fig. 3. What 
is less clear is the onset of this behaviour - there are no large jumps in the estimated yc 
at Up, such as occur for at least some of the approximants to the square lattice series. It 
is fair to say however, that all the plotted approximants (and others we have not plotted) 
do begin to increase around, or just above, the calculated value of Up for the $^ graphs, 
namely 1/9. 

On balance therefore, the analysis of the high-field series for the Ising model on planar 
$^ random graphs also lends support to the existence of a Kertesz line. Since clusters of 
spins behave in the same manner on ensembles of planar random graphs as they do on the 
square lattice (the Wolff |T^ or Swendsen-Wang |T^ algorithms, for instance, still do an 



excellent job of reducing critical slowing down for Ising models on ensembles of random 
graphs) this is reassuring. It would have been interesting, but beyond the scope of the 
resources available to us, to extend the order of the high-field series to see what effect 
this had on the estimates of yc- It would also be an interesting exercise to extend the 
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now vintage square lattice series for further comparison, bearing in mind that the clearest 
signal for a jump in i/c was only seen in the highest available approximants. 
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Figure 1: A schematic drawing of the Kertcsz hne in the M, m = exp(— 4/?) plane. The 
phase co-existence hne is shown in bold, the Kertesz line dotted and the critical Uc and 
percolation Up points are marked. 



9 




Figure 2: The estimated radius of convergence for the high- field series of ^ on the square 
lattice plotted against u. The different symbols represent the different approximants used. 

Note that the jump is seen in both the [7/7] (star) and [6/7] (diamond) approximants, 
whereas the [6/6] (cross) and [5/5] (triangle) approximants give complex roots for some 

u ~ Up. 
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Figure 3: The estimated radius of convergence for the high- field series of ^ on 
graphs plotted against u. The symbols are the estimates obtained from [13/13] (diamond), 
[14/14] (circle) and [15/15] (triangle) Pade approximants, and the [7/7] approximant for 
the square lattice is also shown for comparison using the starred symbol. 
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